Negative thermal expansion and antiferromagnetism in the actinide oxypnictide 

NpFeAsO 
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A neptunium analogue of the LaFeAsO tetragonal layered compound has been synthesized and 
characterized by a variety of experimental techniques. The occurrence of long-range magnetic order 
below a critical temperature Tn = 57 K is suggested by anomalies in the temperature-dependent 
magnetic susceptibility, electrical resistivity, Hall coefficient, and specific heat curves. Below Tn, 
powder neutron diffraction measurements reveal an antiferromagnetic structure of the Np sublattice, 
with an ordered magnetic moment of 1.70 ± 0.07 /j,b aligned along the crystallographic c-axis. No 
magnetic order has been observed on the Fe sublattice, setting an upper limit of about 0.3 \xb for 
the ordered magnetic moment on the iron. High resolution x-ray powder diffraction measurements 
exclude the occurrence of lattice transformations down to 5 K, in sharp contrast to the observation of 
a tetragonal-to-orthorhombic distortion in the rare-earth analogues, which has been associated with 
the stabilization of a spin density wave on the iron sublattice. Instead, a significant expansion of 
the NpFeAsO lattice parameters is observed with decreasing temperature below Tn, corresponding 
to a relative volume change of about 0.2% and to an invar behavior between 5 and 20 K. First- 
principle electronic structure calculations based on the local-spin density plus Coulomb interaction 
and the local density plus Hubbard-I approximations provide results in good agreement with the 
experimental findings. 

PACS numbers: 75.50.Ee, 74.70.Xa, 65.40.Dc, 61.05.fm 



I. INTRODUCTION 

Rarely in physics has a single family of com- 
pounds generated as much interest as the iron-pnictide 
superconductors^— . Research in this field is being 
conducted at such a hectic pace that no fewer than 
2500 articles have cited the paper that reported the 
discovery of superconductivity in LaFcAsOi-^F^. 
The reason for this frenzy of activity is clear: the su- 
perconducting transition temperatures are very high, 
second only to the cuprates. 

The replacement of the rare-earth (R) species in 
the iZFeAsO "1111" phase with Np, without change 
to the room temperature structure, and the ability to 
subsequently replace Np with heavier actinides: Pu 
and Am, represents a real possibility to control cor- 
relations (from weak to strong, Np to Am) and to 
potentially extend the already rich physics of the iron 
pnictides. 

The parent iJFeAsO "1111" compounds have so 



far involved rare-earth 3 + ions in the oxide layers. 
These have ionic radii ~ 1 A and have partially filled 41 
shells, which are atomic-like in nature. The influence 
of the oxide layer, in the behavior of these materials, 
beyond being a 'spacer', is a seldom addressed ques- 
tion. In fact, from a materials point of view this is 
a difficult problem to investigate, since elements that 
possess similar chemical properties to the rare earths, 
both in terms of bonding and size, are hard to find. 
For ideal candidates one needs to descend to the bot- 
tom of the periodic table, specifically, the transuranic 
actinide elements, Np, Pu and Am. These have ionic 
radii of almost identical magnitude to the majority 
of the lanthanides, and are commonly present in the 
trivalent state in similar compounds. Not only do they 
represent an alternative to the lanthanides per se, but 
they can also be used to study the effect of the degree 
of localization of the / electron states (in this case 
5f) on the observed properties. This is a commonly 
adopted strategy in the actinide community when in- 



2 



vcstigating new compounds; the behavior of the ac- 
tinide metals as one traverses the period from ura- 
nium to americium changes from itinerant, transition 
metal-like electronic behavior to localized, atomic-like 
behavior. 

Here, we report on the synthesis and characteri- 
zation of NpFeAsO, an isostructural Np analogue of 
LaFcAsO. Magnetic, electrical transport, and specific 
heat measurements have been used to characterize 
its macroscopic physical properties, whereas neutron 
and x-ray diffraction have been used to monitor the 
evolution of the magnetic and crystallographic struc- 
tures with decreasing temperature. The experiments 
show the development of antifcrromagnetic order on 
the Np sublattice below Tn = 57 K, and the absence 
of observable crystallographic distortions down to 5 
K. This is at variance with the behaviour of the rare- 
earth analogues i?FcAsO, where the occurrence of a 
tetragonal-to-orthorhombic lattice transformation is 
observed well above the critical temperature of the 
rare-earth sublattice. On the other hand, the mag- 
netic transition in NpFeAsO is accompanied by a large 
expansion of the unit cell volume below Tn followed 
by an invar behavior below about 20 K. First-principle 
calculations of the electronic structure provide a Fermi 
surface with almost 2-dimensional character, very 
similar to the one previously reported for i?FeAsO, 
and predict an antiferromagnetic ground state with a 
staggered magnetization along the tetragonal c-axis. 
The anomalous thermal expansion is qualitatively de- 
scribed within a local-density plus Hubbard-I approx- 
imation. 



II. EXPERIMENTAL 

The first and not inconsiderable hurdle one has to 
overcome when studying these materials is sample 
synthesis. This is non-trivial in the case of the ra- 
dioactive and toxic actinides. One point of caution 
at this stage is that we have started our study with 
Np, whereas it is most common to begin with U. For 
the case of the oxypnictides, replacing the rare-earth 
element by U is highly unlikely to be successful, given 
that U can be either tri- or tetravalent, and is the lat- 
ter in the vast majority of oxygen containing U-bascd 
compounds. 

Polycrystalline samples of NpFeAsO were synthe- 
sized at the Institute for Transuranium Elements 
(ITU) by solid state reaction using Fe304 (Alfa Ae- 
sar 99.997%), elemental Fe (Alfa Aesar 99.998%), and 
crystals of NpAs as starting materials. The thor- 
oughly mixed powder was pressed into a pellet, sealed 
in an evacuated silica ampoule and heated at 900 °C 
for 48 hours, before finally being furnace quenched. 
All operations were carried out in a radioprotected 
glovcbox with low oxygen and water concentrations. 

High resolution X-ray diffraction data were mea- 



sured as a function of temperature (5 < T < 300 K) 
on the powder diffraction beamline, ID31 at the Euro- 
pean Synchrotron Radiation Facility (ESRF), France, 
using an incident beam energy of 35 keV. Powderized 
sample for this experiment was doubly encapsulated, 
being put inside a kapton tube (50 /im) and inserted 
in a plexiglass hollow cylinder with a wall thickness of 
200 fj,m. 

The magnetic susceptibility measurements were 
performed at ITU on a Quantum Design MPMS- 
7 SQUID magnetometer, in applied magnetic fields 
[IqHi = 5 T and fioH2 = 7 T, where we define 
XDC = [M(iJ 2 )-M(ifi)]/(tf 2 --ffi). Electrical resis- 
tivity, Hall effect and heat capacity were measured us- 
ing a Quantum Design Physical Properties Measure- 
ment System (PPMS-9). The electrical resistivity was 
determined using a standard 4-probe DC technique, 
with four 0.05 mm diameter platinum wires glued to 
the sample using silver epoxy (Epotek H20E). The 
Hall resistance (Rh) was determined by voltage mea- 
surements Vh under applied magnetic field at +9 and 
-9T. The field response V?/(B) at fixed temperatures 
has been measured to confirm results obtained when 
ramping in temperature. The heat capacity was mea- 
sured using a standard relaxation calorimetry method, 
where the data were corrected for the contribution of 
the Stycast encapsulation material (2850 FT) by using 
an empirical relation determined previously. 

Neutron powder diffraction was performed on the 
high-intensity two-axis D20 diffractometer at the In- 
stitut Laue Langevin (ILL), France, in the high flux 
mode, with an incident wave- length of 2.42 A, using 
a (002) pyrolitic graphite monochromator. About 250 
mg of NpFeAsO powder was used, within a triple- 
walled Al container. Data were collected at 5, 30, 90, 
190, and 290 K using the ILL trans-uranic samples or- 
ange cryostat. Variable temperature ramp data were 
also collected between 30 and 90 K. 



III. RESULTS AND DISCUSSION 

Upon the successful synthesis of NpFeAsO, our ini- 
tial step was to check the crystal structure and sam- 
ple purity using x-ray powder diffraction. At room 
temperature this indicates that this Np compound is 
an isostructural analogue to the celebrated rare-earth 
iron oxypnictide family iZFeAsO, as demonstrated by 
the successful refinement of the high resolution data 
shown in Figure [T] to the ZrCuSiAs-type structure. 
Extra lines appearing in the diffraction pattern cor- 
respond to Np02, which is the only observable impu- 
rity phase (about 8% in weight). On cooling to 5 K 
no orthorhombic distortion was observed within the 
resolution of our data (inset (a)), in contrast to all of 
the i?FeAsO compounds (such as SmFeAsO in inset 
(b)). The absence of the crystallographic distortion 
immediately prompts a question regarding the nature 
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FIG. 1. (Color online) Fit (solid blue line) of the refined structural model (Rietveld method) to the high-resolution room- 
temperature x-ray powder diffraction data (red circles) for NpFeAsO, confirming the tetragonal P4/nmm ZrCuSiAs-type 
structure. The data were collected on the ID3I (ESRF) high resolution x-ray powder diffractometer with a wave-length 
of A = 0.35 A. Tick marks are shown for two refined phases: NpFeAsO (upper, blue) and an 8% Np02 impurity (lower, 
green). Insets (a) and (b) make a comparison of data for NpFeAsO and SmFeAsO^, where only SmFeAsO shows the 
splitting of the room temperature (red) tetragonal (P4/nmm) (111) Bragg reflection into the orthorhombic {Cmma) 
(201) and (021) reflections at 5 K (blue). Panels (c) and (d) show the variation in the a and c lattice parameters, and of 
the unit cell volume as a function of temperature, as extracted from the x-ray powder diffraction measurements, revealing 
a marked negative thermal expansion below T ~ 60 K in NpFeAsO. The red line through the volume data is a fit to the 
second order of the Griineisen approximation for the zero pressure equation of state, as discussed in the text. 



of the magnetism in NpFeAsO, since it is believed that 
in the rare-earth compounds the distortion is intrin- 
sically related to the development of a spin density 
wave on the iron, breaking the rotational symmetry 
within the plane. Inspection of Figure [5] reveals that 
bulk thermodynamic magnetic and transport property 
measurements give no indication of an anomaly within 
the standard temperature range for the i?FeAsO spin 
density wave ordering (120 — 140 K). Instead, anoma- 
lies are only observed at 25 K, corresponding to the 
phase transition in the Np02 impurity phased, and 
significantly at ~ 60 K. 

Even more striking and unusual, however, is the ob- 
served variation of the lattice parameters as a function 
of temperature. As shown in FiguresQJc) an d (d), the 
standard thermal expansion behavior stops abruptly 
below T ~ 60 K, and as the temperature is decreased 
the crystallographic cell expands. The very accurate 
data obtained at the ESRF ID31 high- resolution pow- 
der diffraction beamline allows us to study the effect 
along both principal crystallographic directions in de- 
tail. Although the temperature dependences for both 
a and c lattice parameters, show clear minima, T m i n , 
the c-axis minimum is about 10 K higher than in the 
a-axis direction. Below T m i n the c lattice parameter 
increases monotonically and saturates for T < 2QK. 



On the contrary, the a lattice parameter reaches a 
maximum at around 30 K and then starts to decrease 
with decreasing temperature. The relative lattice size 
change differs by a factor of 4, the one in the c-axis di- 
rection being the greater. The above conclusions are 
not affected by the presence of the Np02 impurity, 
whose lattice parameter shows a regular temperature 
behavior with a contraction at the onset of electric 
quadrupole order (To = 25 K)£. 

Using the second order Griineisen approximation 
for the zero pressure equation of stat o 10 ' 11 , the vol- 
umetric data was modeled according to: 



where U is the internal energy, which we have calcu- 
lated following the Debyc approximation: 

/ T \ 3 r e D /T 3 d 
U(T } = M kB T(§-) I ±*. (2) 

N is the number of atoms in the unit cell, ks is the 
Boltzmann constant, and 8d is the Debye tempera- 
ture, obtaining @ D = 287 K, Q = 0.52 eV, b = 260 



4 




FIG. 2. (Color online)Bulk physical properties of NpFeAsO. (a) The magnetic susceptibility \ = AM/ AH data (closed 
circles) and d(xT)/dT (solid line). The inset shows a Curie- Weiss fit to the inverse susceptibility from which an effective 
paramagnetic moment of 2.78 ± 0.06/ib and a Curie- Weiss temperature of —128 ± 2 K are obtained, (b) The electrical 
resistivity measured in H — T. The inset highlights the magnetically ordered region, showing two fits to the data, 
indicating that the low temperature resistivity does not vary simply as T 2 (red line), but instead an improved fit is 
obtained on including an antiferromagnetic gap (blue line), (c) The Hall coefficient as a function of temperature. The 
inset shows the least squares fit described in the text, (d) The heat capacity, including an inset showing a fit to the low 
temperature data. 



and V = 123.78 A 3 , and hence a negative thermal 
expansion of cj = AV/V = 0.2%. 

The presence and form of the anomaly in the mag- 
netic susceptibility (Fig Ufa)) at ~ 57 K, implies 
the onset of antiferromagnetic ordering, and there- 
fore that the anomalous thermal expansion is a con- 
sequence of a strong coupling between the magnetism 
and the lattice. In order to estimate the precise value 
of the Neel temperature, we followed Ref.—, and plot- 
ted d(xT)/dT (solid line) on the same temperature 
scale. The maximum of d(xT)/dT is observed at 
Tn = 57 K. Above ~ 100K, the susceptibility curve 
has a Curie- Weiss behavior, with an effective param- 
agnetic moment /U e g = 2.78 ± 0.06 /ib and a Curie- 
Weiss temperature Ocw = —128 ± 2 K. The above 
values have been obtained by fitting the inverse sus- 
ceptibility curve to a straight line, taking into account 
the contribution of the Np02 impurity phase by scal- 
ing the data reported in Ref. [H. 

In the intermediate coupling scheme, fi c g of Np 3+ 
is 2.755 fiB, and hence our results imply that the nep- 
tunium is trivalent, as is the case for the rare-earth 



species in the RcFcAsO isostructural compounds. The 
magnetic moment on the iron is likely to be small 
given that the majority of the magnetism can be at- 
tributed to the Np moments. The negative Curie- 
Weiss temperature suggests the presence of antifer- 
romagnetic interactions. An additional anomaly is 
present in the data at T = 25 K. This corresponds 
to the magnetic triakontadipolc phase transition of 
NpOz^— , which is present as an impurity phase in 
our sample, as detected in the x-ray diffraction pat- 
tern. 

Figure Hfb) shows the resistivity, revealing a sharp 
drop at the magnetic transition. The low tempera- 
ture resistivity can be described by the presence of 
AF interactions at high temperatures with the open- 
ing of a magnetic gap near the magnetic transition 
temperature^, as demonstrated by the fit in the in- 
sert to the function (blue solid line): 

p{T) = Po + aT 2 + bT(l + 2T/A) cxp~ A / T , (3) 

where p = 70±l^cm, a = 0.27 ± 0.01 /ificmK" 2 , 
b = 20 ± 1 fiflcmK- 1 and A = 49 ± 1 K. The value 
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of A is comparable to that obtained for NpCoGas 
(A = 55K), which orders antiferromagnetically at 
T N = 47 KiS. 

Whilst the electrical resistivity indicates that 
NpFeAsO is a reasonable metal, the number and na- 
ture of the carriers were further investigated by means 
of the Hall Effect. The temperature dependence of the 
Hall effect of NpFeAsO is shown in Fig.[5Jc). At room 
temperature the Hall coefficient Rh is positive and of 
the order of 2xl0 -9 m 3 C -1 . It increases with de- 
creasing temperature down to the Neel temperature 
Tn = 57 K where Rh(T) exhibits a distinct maxi- 
mum, reminiscent of the maximum in the tempera- 
ture dependence of the magnetic susceptibility (see 
Fig.[2][a)). It has been shown experimentally and the- 
oretically, based on the Anderson periodic model with 
the crystal electric field effect, that in magnetic ma- 
terials the Hall coefficient may be described as a sum 
of two parts: 

R H (T) = R + RsX*(T), (4) 

where i?o is the normal Hall effect due to the Lorenz 
motion of carriers and the second term, the anomalous 
Hall effect, is related to the magnetic scattering^—. 
The reduced susceptibility x*{T) is approximated by 
\(T)/C where C is the Curie- Weiss constant^. In 
the case of NpFeAsO, a least squares fitting of the 
above equation to the experimental data in the tem- 
perature range 60 — 200 K (see the inset to Fig. HJc)) 
resulted in the values i?o = —9.5 x 10 -9 va^C^ 1 and 
R s = 2.1 x 10~ 6 it^C -1 for the normal and anomalous 
Hall coefficients. The anomalous Hall effect, likely 
caused by the magnetism of the Np ions, is domi- 
nant, and therefore a positive Rh is observed. How- 
ever, the negative value of Rq indicates that elec- 
trons are the dominant carriers, most probably due to 
their higher mobility. The single band model provides 
an estimate for the concentration of free electrons to 
be 6.5 x 10 20 cm' 3 , which should be considered as 
the upper limit of the actual carrier concentration in 
NpFeAsO. Interestingly, the values of the ordinary 
Hall coefficient, as well as its sign, and carrier con- 
centration obtained for NpFeAsO are similar to the 
ones derived for similar i?FcAsO systems based on 
the lanthanide elements 2 ^—. Neither resistivity, nor 
Hall effect measurements on NpFeAsO show anoma- 
lies as observed around 150 K in transport properties 
measurements of iJFeAsO (R = La, Ce, Pr, Nd) 2 ^. 

Figure [UJd) presents the heat capacity data mea- 
sured down to 2K. Two anomalies are clearly iden- 
tified, namely a sharp A- type cusp at ^60 K, that 
coincides with the onset of antiferromagnetic order 
in NpFeAsO, and a shoulder at 25K, associated with 
the multipolar phase transition of Np02 impurity 2 ^. 
The insert shows a fit to C P /T = 7 + (3T 2 in the 
3.5 — 7 K temperature range, where 7T is the elec- 
tronic contribution to the heat capacity, and /3T 3 
(/3 = 12-7T 4 N 'ks /50£>) is the acoustic phonon contribu- 



tion in the low temperature limit of the Dcbye model, 
where Qp> is the Debye temperature, N is the number 
of atoms per formula unit, and fcs is the Boltzmann 
constant. The contribution of the Np02 impurity to 
the specific heat data have been subtracted by scaling 
the data reported for Np02 in ref. [26|. 

In the magnetically ordered state the magnetic spe- 
cific heat term (C mag ) should be taken into account. 
In ref. [HI, the formula that describes the specific heat 
of magnons with an energy gap A in their disper- 
sion relation, C mag = aT x l 2 exp (— A/T), was suc- 
cesfully used in order to fit C mag of NpCoGas. Using 
the value of A = 49 K obtained from the resistivity 
fit, and taking a = 3.7 J mo\~ 1 K~ 3 ^ 2 , as reported 
for NpCoGa5, C mag term can be estimated. In the 
3.5 K to 7 K temperature range, C mag changes from 
6^tJmol~ 1 K~ 1 to 9mJmol~ 1 K _1 , which is 5 and 2 
orders of magnitude smaller than the measured value 
of Cp and does not influence the fit shown in the in- 
set of Figure [2^d) , where no deviations from linear 
behavior are observed. The values of 7 and Qd ex- 
tracted from the fit arc 75 ± 1 mJ mol -1 K -2 , and 
178 ± 2 K respectively. The obtained Debye tempera- 
ture is not in good agreement with the value obtained 
from the fit of volumetric data. The observed up- 
turn in Cp I T below T = 3 K is caused by the nuclear 
Schottky anomaly, which is commonly observed in Np 
compound o 18 i 27 . In order to calculate the magnetic 
entropy, the phonon contribution to the specific heat 
must be subtracted. However, we can not use reported 
data for parent i?FeAsO because all of them undergo 
a structural phase transition such that at low temper- 
ature they posses the orthorhombic crystal structure. 
According to our knowledge, despite of a lot of ef- 
fort, synthesis of UFeAsO and ThFeAsO has not been 
successful. Therefore, the entropy removed upon mag- 
netic ordering cannot be calculated at this moment. 

Given both 7 and x(0 K), which was obtained from 
the interpolation of the C - W fit to T = K, the 
Wilson ratio, 

Rw = (^k B ) 2 x(0K)/3 M 2 ff7 , (5) 

can be estimated for NpFeAsO, yielding a value for 
Rw of 1.3. This is slightly higher than the free elec- 
tron value (Rw = 1) an d less than expected for heavy 
fermion compounds (Rw = 2). 

Information about the magnetic structure was ob- 
tained from neutron powder diffraction. On cooling 
below ~ 57 K additional Bragg reflections appeared 
(Figure [3]). Their absence in the high resolution x-ray 
powder diffraction data, and their temperature depen- 
dence reflects their magnetic origin. 

Rietveld modeling of the neutron data were per- 
formed using the GSAS program 2 ^ with EXPGUI 2 ^. 
The magnetic form factor assuming spherical symme- 
try (jo and j'2) and the bound coherent scattering 
length of 237 Np 3+ (10.55 fm) were taken from the in- 
ternational tables^. A Landc splitting factor (g) of 
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FIG. 3. (Color online) Variable temperature neutron 
diffraction surface across Tn, showing the appearance of 
the structurally forbidden (100) reflection and the increase 
in intensity of the (101) and (111) reflections due to mag- 
netic ordering. 

0.6 was used. Several parts of the diffraction patterns 
were excluded, due to the aluminium container and 
other encapsulation and sample environment artifacts. 

The 5 K and 30 K data sets were fitted with the 
AFM PA/n'm'm' model, whilst the remaining sets 
were fitted with a paramagnetic P4/nmm model, the 
results of which are shown in Table I. To confirm that 
all the magnetic Bragg intensities were being mod- 
eled correctly, a purely magnetic fit was carried out 
to the difference between the 30 K and 90 K data sets 
(Figured]). A list of reflections, with the measured 
and calculated structure factors are shown in Table 
HI1 These peaks indexed using the PA/n'm'm! mag- 
netic Space Group, correspond to the Np moments 
being aligned along the c-axis, with a ferromagnetic 
coupling within the basal plane, and an antiferromag- 
netic coupling between the planes (see the inset to 
Figure |4]) , reminiscent of the low temperature order- 
ing of the Sm moments in SmFeAsO 31 ' 32 . 

No Fc moment is observed within the limits of the 
experiment. At 5 K, the ordered magnetic moment of 
the Np ions is of 1.70±0.07^tB, in excellent agreement 
with Mdssbauer spectroscopy estimations^ 3 .. 

It is worth noting, that for the parent i?FcAsO 
"1111" compounds, the magnetic rare-earth moments 
order with the moments varying from 0.83 to 1.55 for 
CeFeAsO and NdFcAsO, respectively^. The high- 
est Neel temperature, of the rare-earth sublattice, in 
the i?FeAsO "1111" series, T N (Pr)=14 K observed for 
PrFcAsO 3 ^, is four times lower than Tn = 57K esti- 
mated for NpFeAsO. 

To explore further the electronic and magnetic char- 
acter of NpFeAsO we have performed first-principle 
local spin density approximation (LSDA) and LSDA 
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FIG. 4. (Color online) Plot of the Rietveld fit to the dif- 
ference of the 90 K and 30 K neutron powder diffraction 
data sets using a P4/n'm'm' magnetic model. Portions of 
the difference spectra where strong nuclear reflections led, 
due to thermal expansion, to strong up-down features were 
excluded from the fit. The right inset shows the constant 
presence of background features at all temperatures in the 
mid momentum transfer region. The left inset shows the 
proposed antiferromagnetic structure. 



plus Coulomb-U (LSDA+U) calculations. In all cal- 
culations we used the crystal structure parameters 
determined experimentally at room temperature, as- 
suming the tetragonal ZrCuSiAs structure. Gaining 
inspiration from the neutron scattering data, we as- 
sumed that the magnetic and crystallographic unit 
cells coincide. Furthermore, we assumed that the 
magnetic moment is aligned along c-[001]-axis, and 
considered non-magnetic (NM), ferromagnetic (FM), 
and anti-ferromagnetic (AF) arrangements for the Np 
moments. There were no magnetic moments set ini- 
tially on the Fe atoms. In all calculations we assumed 
stoichiometric NpFeAsO and did not consider any de- 
fects or impurity Np02 phase. The effect of defects 
and impurities is left for further consideration. 

We used an in-house implementation of the full- 
potential linearized augmented plane wave (FP- 
LAPW) method 3 ^. This FP-LAPW version includes 
all relativistic effects (scalar-relativistic and spin-orbit 
coupling), and relativistic implementation of the ro- 
tationally invariant LSDA+U 3 ^. In the FP-LAPW 
calculations we set the radii of the atomic spheres to 
2.75 a.u. (Np), 2.2 a.u. (Fe, As), and 1.6 a.u. (O). 
The parameter Rn p x K max = 9.625 determined the 
basis set size, and the Brillouin zone (BZ) sampling 
was performed with 405 k points. For the neptunium 
/ shell, Slater integrals of F a = 3.00 eV, F 2 = 7 A3 eV, 
F 4 = 4.83 eV and F 6 = 3.53 eV were selected to spec- 
ify the Coulomb interaction 3 ^. They correspond to 
commonly accepted values for Coulomb U — 3 eV 
and exchange J = 0.61 eV parameters. The depen- 
dence of Np-atom mj on the choice of Coulomb-?/ in 
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TABLE I. Crystallographic parameters as a function of temperature derived from Rietveld refinement of the structural 
model to match the neutron powder diffraction data from D20. ZrCuSiAs structure, space group P4/nmm (origin choice 
2) with Np, Fe, As, and O on the 2c, 2b, 2c, and 2a Wyckoff sites, respectively. The Np magnetic moment (Np m z ) 
refined in the PA/n'm'm' magnetic space group in the low temperature phase is also presented here. 



T 1 K 




5 


30 


90 


190 


290 


a/ A 




3.85525(15) 3.85512(15) 3.85530(15) 3.85851(16) 3.86235(17) 


c / A 




8.3374(6) 


8.3371(6) 


8.3324(6) 


8.3466(7) 


8.3663(7) 


Np z 




0.1507(6) 


0.1508(6) 


0.1507(6) 


0.01512(7) 


0.1519(8) 


As z 




0.6742(11) 


0.6742(10) 


0.6755(10) 


0.6749(10) 


0.6738(12) 


Np U iso / A 2 


0.028(2) 


0.029(2) 


0.027(2) 


0.031(3) 


0.036(3) 


Fe U iso 


/A 2 


0.032(3) 


0.032(3) 


0.033(3) 


0.034(3) 


0.034(3) 


As U isc 


/A 2 


0.032(4) 


0.032(4) 


0.035(4) 


0.037(4) 


0.038(4) 


U iso 


/A 2 


0.028(5) 


0.029(5) 


0.027(4) 


0.033(5) 


0.039(5) 


Np m z 


/ 


1.70(7) 


1.66(7) 








-/?wp 




0.0212 


0.0209 


0.0205 


0.0201 


0.0198 


Rp 




0.0155 


0.0152 


0.0151 


0.0146 


0.0141 



a range from 1 eV to 3 eV was checked, and only small 
changes in the rrij moment value were found. 

In LSDA, we found that spin-polarization decreases 
the total energy with respect to the non-magnetic so- 
lution by 0.351 eV/ per formula unit (f.u.) for the FM- 
solution, and by 0.335 cV/f.u. for the AF-solution, 
suggesting an FM-ordered ground state. However, in- 
cluding the Coulomb- U and exchange- J in the fully- 
localized-limit (FLL) LSDA+U method, we obtain a 
total energy for the FM solution which is higher by 
12 meV/f.u. than in the AF solution; on the con- 
trary, if the double counting correction is made in 
the around- mean- field- limit (AMF), the AF solution 
becomes lower in energy than the FM one by 61.2 
meV/f.u. Thus, both flavors of LSDA+U yield an 
AF-ordered ground state, in agreement with the ex- 
periment. 

The spin ms, orbital tul, and total mj magnetic 
moments for the AF calculations with LSDA and 
LSDA+U are shown in Tabic IHII The staggered local 
magnetic moments, which are due to magnetic polar- 
ization of the Np /-shell are formed at the two Np 
atoms. Comparison with the experimentally deter- 
mined moment of 1.7 fiB shows that LSDA is failing 
completely to reproduce the results of neutron scat- 
tering experiments, whilst the AMF-LSDA+U calcu- 
lations overestimate the mj moment value. On the 
other hand, the FLL-LSDA+U calculations yield an 
mj value in reasonable agreement with the experi- 
mental data, and the ratio of the orbital to total mo- 
ment gives a g-value of 0.58, consistent with that for 
a free Np 3+ ion, and the value of 0.6 used in the NPD 
refinements. 

It should be noted that when the total energy dif- 
ference between the FM and the AF solutions is used 
to make a "naive" molecular-field theory estimate of 
the Nccl temperature T N , the FLL-LSDA+U yields 
Tn ~ 47 K in reasonable agreement with the exper- 



imental data. Thus we assume that the LSDA+U 
model with the FLL-double-counting choice gives the 
most appropriate description of the NpFeAsO. 

Although no magnetic polarization on other atoms 
in the unit cell was initially assumed, the presence of 
the exchange splitting at the Np atoms leads to in- 
duced staggered spin moments on the Fe atoms in a 
checker-board AF arrangement. These are fairly neg- 
ligible in LSDA, but reach ~ 0.1-0.25 {1b values in 
LSDA+U. Such a moment on the iron was not visi- 
ble in our neutron powder diffraction measurements, 
but can not be excluded due to the background aris- 
ing from the sample encapsulation and to the high 
absorption cross-section of Np. 

The total and partial (atom, spin and orbital- 



TABLE II. Observed, F 2 , and calculated, F 2 , structure 
factors for the magnetic Bragg reflections in the 30 — 90 
K refinement. Structure factors have been corrected for 
scale and extinction. 



hkl Q (A' 1 ) Fg F 2 



100 


1 


.63 


0, 


9312 





,9128 


101 


1 


.80 


0, 


2374 





,2480 


102^ 


2 


.22 









04335 


111 


2 


,43 


0, 


,4749 





,4442 


112^ 


2 


.75 









,4174 


103^, 


2 


,79 


0, 


09241 





,2043 


113^ 


3 


,23 









02390 


201 


3 


,35 


0. 


3697 





,3177 


104^ 


3 


,43 









06999 


202 


3 


,59 





3496 





3294 


210 


3 


.64 


0. 


3925 





,4317 


211 b 


3 


.72 





1687 





,1349 



a Excluded region. 

b Partially excluded region. 
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LSDA AMF-LSDA+U FLL-LSDA+U NPD 

Np Fe As O Np Fe As O Np Fe As O Np 

m s 2.91 0.00 0.03 0.00 1.53 0.25 0.01 0.00 2.98 0.10 0.02 0.00 2.26 

m L -3.17 0.00 0.00 0.00 -4.23 0.02 0.00 0.00 -5.03 0.01 0.00 0.00 -3.96 

mj -0.25 0.00 0.03 0.00 -2.70 0.27 0.01 0.00 -2.05 0.11 0.02 0.00 -1.7 

TABLE III. Spin ms, orbital rriL, and total mj magnetic moments (/xb), for one (of two) Np, Fe, As, and O atoms in 
anti-ferromagnetic NpFeAsO resulting from calculations with LSDA and LSDA+U=3 eV with AMF- and FLL-double- 
counting flavors. The last column (NPD) provides the values obtained from the neutron powder diffraction studies 
making use of the g-value of for a free Np 3+ ion. 



— Total 

— Np-f 

— Fe-d 
As-p 
O-p 


k 




1 




if 




1 













-7 -6 -5 -4-3-2-10 1 2 3 

Energy (eV) 

FIG. 5. Electronic structure calculations for NpFeAsO: 
the relativistic spin- and orbital-resolved DOS (per unit 
cell) including the partial DOS 



resolved) density of states (DOS) are shown in Fig. 5 
calculated with FLL-LSDA+U= 3 eV. The DOS near 
the Fermi energy (Ep) has mostly Fe-c? character, 
while the As-p and O-p states are mostly located at 
2-7 eV energy interval below Ep. The Np-/ states 
are split by the exchange interaction. The resulting 
DOS at Ep of 4.8 states/eV corresponds to the non- 
interacting value of the Sommcrfcld coefficient 7= 5.6 
mJ mol _1 K -2 which is substantially lower than the 
experimental value of 75 mJ mol _1 K -2 . However, 
the theoretical value of the Sommerfeld coefficient is 
expected to increase due to electron mass enhance- 
ment caused by dynamical electron interactions and 
electron-phonon coupling. 

The band structure and the Fermi surface (FS) 
are shown in Fig. 6. The FS consists of five sheets, 
each of them doubly degenerate. Examination of the 
band structure shows that FS-1-3 sheets are hole-like, 
and centered at the T-point. The FS-4 and FS-5 are 
electron-like and centered at the M-A line. Note the 
fairly two dimensional character of the FS, and the 
strong resemblance to those previously presented for 
fiFcAs O 39 ' 40 . It is seen that most of the states in the 
vicinity of the FS are located near T-[0, 0,0] and M- 
[tt, tt, 0] k-points in the BZ suggesting the possibility 



of s± superconducting pairing mechanisms^. 

A striking feature of NpFeAsO which makes it dif- 
ferent from the rare-earth-based counterparts, is the 
absence of the orthorhombic structural distortion as- 
sociated with the magnetic ordering. In order to un- 
derstand this, we performed magnetic anisotropy en- 
ergy (MAE) calculations, rotating the staggered AF 
magnetization from the c-axis to the a-axis direction. 
What we found is that the total energy E to t differ- 
ence E tot (a) — E to t{c) between these two directions 
of the magnetization is of « 30.0 meV/f.u. This 
means that there is a strong positive uniaxial MAE in 
NpFeAsO. This MAE keeps the staggered AF magne- 
tization along the tetragonal c-axis and assists in the 
prevention of any distortion in the a — 6-plane. 

Finally we turn to the discussion of the negative 
thermal expansion (NTE). The most recent theory of 
the Invar effect^ is based on the use of a disordered 
local moment (DLM) approximation. The paramag- 
netic state in the DLM approximation is treated as 
a disordered pseudo-alloy with equal concentration of 
randomly oriented "up" and "down" local moments. 
It is usually implemented by making use of alloys the- 
ory in a coherent-potential approximation. Regret- 
tably, it is not currently possible to use the DLM ap- 
proximation together with the FP-LAPW basis used 
in this work. 



Method LDA LDA+HIA FLL-LSDA+U 

V eq , (A 3 ) 111.50 116.11 119.92 
B, (MBar) 1.92 1.86 1.79 

TABLE IV. The equilibrium volume V eq (A 3 ) for 
NpFeAsO resulting from paramagnetic LDA+HIA and 
anti-ferromagnetic LSDA+U calculations. 

Therefore, here the LDA -(-Hubbard I (HIA) approx- 
imation, as implemented in Reff43l is used to treat 
the Np-atom local moment paramagnetic phase. The 
same values for the Slater integrals for Np atom as in 
LSDA+U were used in LDA+HIA calculations (which 
correspond to a Coulomb U = 3 eV and an exchange 
J = 0.61 cV), and the inverse temperature /3=40 
eV -1 was chosen. The use of atomic-like LDA+HIA 
is justified for NpFeAsO at high temperature since 
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FS 1 FS-2 FS-3 FS-4 _ FS-5 




FIG. 6. (Left) Band structure for NpFcAsO;(right) Fermi surface for NpFeAsO, calculated making use of the 
FLL+LSDA+U=3 eV calculations. 



the experiments presented above show a well local- 
ized character of the Np-atom f-shell above the Neel 
temperature. Furthermore, the effective paramag- 
netic moment of 2.84 [Lb per Np atom calculated in 
LDA+HIA, is in very good agreement with the mo- 
ment of 2.78±0.06 fiB obtained from susceptibility 
measurements. 

The total energy versus volume was calculated for 
LSDA+U and LDA+HIA, and the equilibrium vol- 
umes {V eq ) were obtained using a 3rd order polyno- 
mial fit. The c/a-ratio and the internal positions of 
Np and As atoms were kept fixed in these calculations. 
The resulting V eq for LSDA+U is reasonably close 
(with 3.3 % difference) to the experimental V eq = 124 
A 3 at T= K, confirming an adequate choice of the 
computational model for NpFeAsO. However, the cal- 
culated w of 3.1%, whilst agreeing in sign, is about ten 
times bigger than the experimental value (u> = 0.2%). 

This quantitative disagreement between theory and 
experiment is traced to the incorrect description of 
spin fluctuations on the Fe atoms in LDA+HIA cal- 
culations. The Hubbard-I correction is included for 
the Np atom only, and describes the local moment 
paramagnetic phase in Np-sublattices. The Fe-atom 
sublattice is treated within LDA, with no spin polar- 
ization (zero magnetic moment). This approximation 
is not sufficient for the Fe-atom sublattice in the para- 
magnetic phase of NpFeAsO. Note that it bears close 
similarities to what has been found in Ref. 13, for the 
ui calculations in ferromagnetic RE(=Gd,Dy,Er)Co2 
compounds. The FM-state was treated with the so- 
called "open-core" approximation for RE-atom, and 
the DLM approximation was used for PM-phase. The 
calculated u> was found to be about five times bigger 
than in the experiment, and no spin polarization on 
the Co-atoms was found in the DLM calculations. By 
forcing a FM-moment on the Co-atom sublattice (ap- 
plying DLM to the RE-sublattice only), the results 
improved quantitatively, and the value of oj reduced 



substantially. 



IV. CONCLUSIONS 

In conclusion, we successfully synthesized 
NpFeAsO, one of the few actinide-based iron- 
oxypnictide compound (the other reported oxypnic- 
tides containing actinide elements are U2Cu2As30^ 
, UCuPOi&il, and NpCuPOil). Although one may 
expect the physical properties of NpFeAsO to be 
similar to the lanthanide 1111 analogues, in reality 
this system behaves differently. In particular, we 
do not observe any structural transition, which is 
in agreement with the fact that neither transport 
properties (resisitivity and Hall effect) nor specific 
heat measurements reveal high temperature anoma- 
lies associated with an in-plane magnetic ordering 
in the Fe-As layer. We suggest that the lack of the 
orthorhombic structural distortion is caused by a 
strong positive uniaxial magnetic anisotropy energy. 

Our results reveal that NpFeAsO exhibits an 
anomalous thermal expansion which is intimately 
linked to a giant magneto-expansion occurring as the 
system undergoes antifcrromagnctic ordering. Neg- 
ative thermal expansion was observed in PrFeAsO^, 
and we believe that it is a common feature in the 1111 
family with a magnetic trivalent ion. However the ef- 
fect in NpFeAsO is at least 20 times stronger and is 
clearly associated with the antiferromagnetism of the 
Np sublattice. 

This invar effect is usually observed in transition 
metal disordered alloys (e.g. Fe-Ni,Fc-Pt42), and 
rare-earth ferromagncts (e.g. GdCu^) and anti- 
ferromagnets (e.g. PrFeAsO^). It was also found 
in a-U metal^ and later this effect was explained by 
the onset of a charge density wave. Even more strik- 
ing is the observation of a NTE in 6- Pu£2,, where the 
absence of magnetic ordering necessitates an alterna- 
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tive mechanism to that based on two spin states with 
different volumes, as for invar alloy a 53 ' 54 . Instead in 
that case it is proposed that the negative thermal ex- 
pansion is a signature of criticality arising from prox- 
imity to a localization-delocalization transition^. Al- 
though uranium and plutonium have rather overshad- 
owed their neighbor, neptunium^, here, in reporting 
the observation of NTE in NpFcAsO, neptunium will 
finally take centre-stage. To the best of our knowl- 
edge, NpFeAsO is the first actinide-based antiferro- 
magnetic material, in which negative thermal expan- 
sion has been observed. This behavior originates from 
the magnetic order, but more exotic explanations, e.g. 
possible criticality at the interface between localised 
and itinerant electrons, should also be taken into con- 
sideration. 
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